Robust Regularity in 7-Soft Nuclei and its Microscopic Realization 
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7-softness in atomic nuclei is investigated in the framework of energy density functionals. By 
mapping constrained microscopic energy surfaces for a set of representative non-axial medium- 
heavy and heavy nuclei to a Hamiltonian of the proton-neutron interacting boson model (IBM- 
2) containing up to three-body interactions, low-lying collective spectra and transition rates are 
calculated. Observables are analyzed that distinguish between the two limiting geometrical pictures 
of non-axial nuclei: the rigid-triaxial rotor and the 7-unstable rotor. It is shown that neither of 
these pictures is realized in actual nuclei, and that a microscopic description leads to results that 
are almost exactly in between the two geometrical limits. This finding points to the optimal choice 
of the IBM Hamiltonian for 7-soft nuclei. 

PACS numbers: 21.10.Re,21.60.Ev,21.60.Fw,21.60.Jz 



Like many other quantum systems, atomic nuclei dis- 
play a variety of geometrical shapes that reflect deforma- 
tions of the nuclear surface arising from collective motion 
of many nucleons [l[ . Shapes of most non-spherical nu- 
clei are characterized by axially-symmetric quadrupolc 
deformations - prolate or oblate ellipsoids. There are, 
however, many nuclei in which axial symmetry, i.e., the 
invariance under rotation around the symmetry axis of 
the intrinsic state, is broken. The precise description 
of axially asymmetric shapes and the resulting triax- 
ial quantum many-body rotors remains open questions 
in nuclear physics and, since they are also being devel- 
oped for other finite quantum systems like polyatomic 
molecules Q, presents a topic of broad interest. 

Quadrupolc shape deformations can be described in 
terms of the polar deformation parameters /? and 7 
The parameter j3 is proportional to the intrinsic 
quadrupole moment, and the angular variable 7 speci- 
fies the type of the shape. The limit 7 = corresponds 
to axial prolate shapes, whereas the shape is oblate for 
7 = 7r/3. Triaxial shapes are associated with interme- 
diate values < 7 < ir/3. The latter have been in- 
vestigated extensively using theoretical approaches that 
are essentially based on the rigid-triaxial rotor model 
of Davydov and Filippov |3|, and the 7-unstable rotor 
model of Wilets and Jean [4| . The former assumes that 
the collective potential has a stable minimum at a par- 
ticular value of 7, whereas in the latter the potential is 
independent of 7 and thus the corresponding collective 
wave functions are extended in the 7 direction. 

However, presumably all known axially-asymmetric 
nuclei exhibit features that are almost exactly in be- 
tween these two geometrical limits, characterized by the 
energy-level pattern of quasi-7 band: relative locations 
of the odd-spin to the even-spin levels. As the two mod- 
els originate from different physical pictures, the ques- 
tion of whether axially-asymmetric nuclei are 7 rigid or 
unstable has attracted considerable theoretical interest 



[H • The present Letter addresses this question from 
a microscopic perspective, and identifies the appropriate 
Hamiltonian of the interacting boson model (IBM) Q 
for 7-soft nuclei, consistent with the microscopic picture. 
We thereby provide a solution to the problem concerning 
the energy-level pattern of the odd-spin states. 

At present the most complete microscopic description 
of ground-state properties and collective excitations over 
the whole chart of nuclides is provided by the frame- 
work of energ y d ensity functionals (EPFs). Both non- 
relativistic and relativistic [H, [H| EDFs have 

successfully been employed in numerous studies of shape 
phenomena and the resulting complex excitation spectra 
and decay patterns fl4H17l |. The starting point is usu- 
ally a constrained self-consistent mean-field calculation 
of the energy surface with the mass quadrupole moments 
as constrained quantities @. This is illustrated in the 
first row of Fig. [TJ where we display the self-consistent 
quadrupole energy surfaces of 134 Ba (a) and 190 Os (b) in 
the (3 — 7 plane. The constrained energy surface of 134 Ba 
is calculated using the relativistic Hartrec-Bogoliubov 
model [HI with the DD-PC1 [lj] functional, and that of 
190 Os employing the Hartree-Fock plus BCS model [lj| 
with the Skyrme functional SkM* [201 . These function- 
als are representative of the two classes - relativistic and 
non-relativistic EDFs, and will be used throughout this 
work to demonstrate that the principal conclusions do 
not depend on the particular choice of the EDF. One no- 
tices in Figs.fTJa) andfljb) that in both cases the energy 
surface is very soft in 7, with 134 Ba displaying a nearly 
7-independent picture, whereas a more pronounced rigid 
triaxial shape is predicted for 190 Os with the minimum 
at 7 ps 30°. 

To calculate excitation spectra and transition rates, it 
is necessary to go beyond the mean-field solution and 
explicitly take into account collective correlations. Here 
we employ the IBM to analyze spectroscopic properties 
of 7-soft nuclei. The 0(6) dynamical symmetry Q of 



2 



|(a) 134 Ba 



DD-PC1 




4 d p 



d\s p 



FIG. 1. (Color online) Self-consistent energy surfaces of 134 Ba 
(a) and 190 Os (b) up to 5 MeV in energy from the minima 
with the DD-PC1 and with the Skyrme SkM* functional, 
respectively. The corresponding mapped energy surfaces of 
the IBM with (middle row) and without (lower row) the three- 
body term of Eq. ((2| are plotted. 



IBM presents a good approximation to a system with 7- 
independent collective potential. The geometrical picture 
of the 0(6) limit of the IBM emerges in the coherent-state 
framework J21( , being consistent with the model of Wilcts 
and Jean [4J. The coherent state represents the intrin- 
sic wave function of the boson system, and 0(6) states 
in the laboratory system can be generated by angular 
momentum projection 21|. The triaxial-rotor features 
of the IBM were emphasized already in [22|, [23| , leading 
to the "equivalence" ansatz of the 7-rigid and the 0(6) 
descriptions of the low- lying spectra [24| . 



The present study uses the proton-neutron IBM (IBM- 
2), which includes proton (neutron) monopolc s T (s„) 
and quadrupole (d u ) bosons, representing J T = + 
and 2 + collective pairs of valence protons (neutrons), re- 
spectively [25j]. The number N w (N v ) of proton (neu- 
tron) bosons equals the number of valence proton (neu- 
tron) pairs (particles or holes) , with respect to the nearest 
proton (neutron) closed shell [25|. The following IBM-2 
Hamiltonian is employed: 

#IBM = e(rc<fcr + Tldu) + kQ„ ■ Q v + i? 3 B (1) 

with the G?-boson number operator n^p = <v p ■ d p (p — 



tt,v), and the quadrupole operator Q p 

Xp[d p d p ]( 2 \ The third term H^b on the right-hand side 
of Eq. fl} denotes the three-body boson interaction: 
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The three-body term was introduced and analyzed in the 
IBM-1 framework (without distinction between proton 
and neutron bosons) [H, [27|, but is used here for the 
first time in the microscopic IBM-2 model. In the IBM- 
2 there could be three-body terms with combinations of 
proton and neutron d-boson operators different from the 
one used in Eq. ([2]). However, since the proton- neutron 
quadrupole interaction dominates over the proton-proton 
and neutron-neutron ones for medium-heavy and heavy 
deformed nuclei, the term ^ represents the dominant 
contribution of three-body boson interactions. For each 
p and p' , there are five linearly independent combinations 
in Eq. determined by the value of L = 0,2,3,4,6 
[26l | . However, only the term with L = 3 can give rise 
to a stable triaxial minimum at 7 w 30° [27|] , because its 
expectation value in the classical limit is proportional to 
cos 2 37. We thus consider only the L = 3 in Eq. © and, 
in addition, assume 6>J = #3 = #3. 

The parameters e, k, Yjr V and #3 are adjusted following 
the procedure of Ref. [14| : the microscopic quadrupole 
energy surface, obtained from a mean-field calculation us- 
ing a given EDF, is mapped onto the corresponding boson 
energy surface, i.e., expectation value of Aibm m the co- 
herent state (cf . [28| for details) . The deduced value 
of 03 > varies gradually with boson number: |#3/k| s» 1 
for 1 < N n + N v < 5 and w 0.5 for 5 < N„ + N u < 10. 

Without three-body boson terms the energy expec- 
tation value either has a minimum at 7 = 0° (prolate 
shapes) or 60° (oblate shapes), or is independent of 7 in 
the 0(6) limit. Triaxial minima are obtained only after 
the inclusion of the three-body interaction H^b- This is 
nicely illustrated in Fig.[TJ where the mapped energy sur- 
faces of the IBM are plotted in the middle row (for the full 
IBM Hamiltonian Eq. ((T|) that contains the three-body 
term), and in the lower row (for the IBM Hamiltonian 
without the three-body term). For 190 Os the Hartree- 
Fock plus BCS model with the Skyrme functional SkM* 
predicts a minimum at 7 m 30°, which can only be re- 
produced on the mapped surface corresponding to the 
expectation value of the full IBM Hamiltonian contain- 
ing the three-body term (Fig.fljd)). The contribution of 
this term to the mapped energy surface is in general less 
important when the number of active bosons becomes 
relatively small. Thus for 134 Ba nucleus in Fig. [l|c) the 
minimum is still on prolate axis even when the three- 
body term is included. The IBM Hamiltonian with up 
to two-body terms yields an energy surface that is soft in 
the 7 degree of freedom (cf. Fig. Hff )) , but the minimum 
is on the 7 = 0° axis. We note that while the angu- 
lar variables 7 of the boson energy surface and the con- 
strained microscopic energy surface are identical to each 
other, the axial deformation parameters (3 are related by 
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a constant of proportionality determined by equating the 
corresponding intrinsic quadrupole moments [l4j|. The 
geometrical variable j3 is obtained by multiplying the bo- 
son axial deformation by factors « 0.15 and 0.2 for 134 Ba 
and 190 Os, respectively. 

A distinction between 7-unstable and rigid-triaxial nu- 
clei arises when considering the ratio of excitation en- 
ergies 0: S(J, J - 1, J - 2) = \{E{J) - E(J - 1)} - 
{E(J-l)-E(J-2)}]/E(2i) for the quasi-7 (K n =2+) 
band J 71 ' =2+, 3+, 4+ . . . , etc. The excitation energies 
E(J) arc obtained by diagonalization of the Hamiltonian 
-Hibm, and the quadrupole operators Q p are used in the 
calculation of E2 transition rates, with identical proton 
and neutron boson effective charges. 

For a characteristic set of non-axial medium-heavy and 
heavy nuclei, in Fig. [5] we plot the energy ratios 5(4, 3, 2) 
(a) and 5(5, 4, 3) (b), as functions of the product of pro- 
ton and neutron boson numbers: N n N u . The latter 
quantity reflects the amount of valence proton-neutron 
correlations, and hence the increase of N^N U corresponds 
to an enhancement of collectivity Q ■ In this work we con- 
sider non-axial nuclei in the mass regions A ~ 110, 130 
and 190, whose spectra display signatures of 7-softness. 
The set of nuclei shown in Fig. [2] has been selected so 
that the corresponding values of N^N^ evenly span the 
widest possible range. The IBM excitation spectra have 
been calculated starting from self-consistent mean-field 
energy surfaces that correspond to the two functionals, 
Skyrme SkM* and the relativistic DD-PC1. The two en- 
ergy ratios, calculated with and without the three-body 
term of Eq. © in the IBM Hamiltonian, arc plotted in 
comparison to data [29j , and the predictions of the rigid- 
triaxial rotor model of Davydov and Filippov Q and the 
7-unstable rotor model of Wilets and Jean [|[ . One no- 
tices that for all considered nuclei data can only be re- 
produced with the IBM Hamiltonian that includes the 
three-body term Eq. Both the empirical and calcu- 
lated ratios fall almost exactly in between the limits of 
the 7-unstablc rotor and the rigid-triaxial rotor models: 
the Wilets- Jean limit is -2.00 and the Davydov-Filippov 
limit is 1.67 for 5(4,3,2); the Wilets- Jean model pre- 
dicts 2.50, and Davydov-Filippov -2.33 for 5(5, 4, 3). The 
IBM Hamiltonian with up to two-body terms cannot re- 
produce the empirical values and, in both cases, yields 
energy ratios that are close to the predictions of the 7- 
unstable rotor model. 

While the energy ratios are largely independent of the 
product of boson numbers, the B(E2) systematics re- 
flects the evolution of collectivity. For instance, the ratio 
B(E2;3^ 2+)/B(E2;2+ -> 0+), plotted in Fig.^c), 
gradually increases with N^Ny. For nuclei with typi- 
cally low N„N V « 10), like 132 < 134 Ba and 194 < 196 Pt, the 
7 value on average is close to 0° or 60°. In this case 
the ratio B(E2;3t -> 2+)/S(E2;2+ ->■ 0+) is closer to 
the Wilets- Jean limit (0(6) in the IBM representation) 
of 1.19. As the collectivity evolves with N^N^ 12, 
this B(E2) ratio, calculated with the full IBM Hamil- 
tonian that includes the three-body term, saturates be- 
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FIG. 2. (Color online) The energy ratios (a) 5(4, 3, 2) and (b) 
5(5, 4, 3), and (c) the S(E2; 3+ -> 2+)/B(E2; 2+ -> 0+ ) ratio, 
as functions of the product N^N^, for a characteristic set of 
non-axial medium and heavy nuclei. IBM(3B) and IBM(2B) 
denote results obtained with the IBM Hamiltonians with up to 
three- and two-body terms, respectively. The Skyrme SkM* 
and relativistic DD-PC1 functionals are used. Data are from 
Refs. [29I. [3C|| . and D-F and W-J denote the limits of the rigid- 
triaxial and the 7-unstable (or 0(6)) models, respectively. 



tween the 7-rigid limit of 1.78 and the 7-unstable limit 
of 1.19, in agreement with behavior of the energy ratios 
5( J, J - 1, J - 2). The B(E2) ratio calculated with the 
IBM Hamiltonian with up to two-body terms remains 
close to the 0(6) limit even for large values of N^N^. 

Although the ratios shown in Fig. [2] are calculated us- 
ing two completely different microscopic density func- 
tionals, it appears that the basic features of this analysis 
are not sensitive to the particular choice of the underlying 
EDF. 

In the IBM picture, the number of proton (neutron) 
bosons equals half the number of the corresponding va- 
lence particles or holes (25|. Among the nuclei discussed 
in this Letter, those with relatively large N^N U (^ 12), 
in many of which both N n and N v correspond to hole 
configurations, are more likely to exhibit pronounced 7 
rigidity, compared to systems with low N^N^ (^ 10). In 
most of the latter cases N v {N v ) corresponds to particle 
(hole) configuration, or vice versa. 

The discussion so far has focused on the systematics of 
energy ratios and transition rates. The model, however, 
provides an equally accurate and complete description of 
low-energy excitation spectra in individual nuclei. This is 
highlighted by the level scheme of 190 Os in Fig. [3] Again 
we compare results obtained with IBM Hamiltonians con- 
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FIG. 3. (Color online) Low-lying spectra and 73(E2) values 
(in Weisskopf units) of 190 Os. The bands calculated with 
IBM Hamiltonians with (IBM(3B)) and without (IBM(2B)) 
the three-body term of Eq. ([2} are compared to experimental 
data [1H. The Skyrme SkM* EDF is used, and the boson ef- 
fective charge is adjusted to reproduce the experimental value 
of B(E2;2+ -> 0+). 



taining up to two- and three-body terms to available data 
[29l l3lj |. The full IBM Hamiltonian Hibm reproduces 
both the excitation energies and transition rates for the 
ground-state band and the band built on the state 
(quasi-7 band). We notice the marked effect of the three 
body term on the quasi-7 band: all states are lowered in 
energy but, in particular, the pronounced lowering of the 
odd-spin states, e.g., 3^ and 5+ by 473 keV and 663 keV, 
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respectively, breaks the quasi-degeneracy of the doublets 
(3^~,4j), (5^",6j), etc [33]. These doublets (r-multiplets) 
are characteristic of the 7-unstable 0(6) limit of IBM Q . 
We emphasize that there are no additional adjustable pa- 
rameters in the calculation of excitation energies, that is, 
the parameters are completely determined by the choice 
of the microscopic functional and the mapping procedure. 
Results of similar level of agreement with experiment are 
also obtained in the calculation of spectra of other nuclei 
considered in this study. 

In conclusion, we have investigated the emergence of 7 
softness in atomic nuclei starting from the microscopic 
framework of energy density functionals. For a wide 
range of relevant nuclei certain observables allow us, 
in comparison to microscopic calculations, to differen- 
tiate two limiting geometrical pictures: the rigid-triaxial 
and the 7-unstable rotors. The present analysis clearly 
demonstrates that neither of these pictures is realized 
in actual nuclei. Typical non-axial medium-heavy and 
heavy nuclei lie almost exactly in the middle between 
the two geometrical limits, as a robust regularity. In 
the IBM framework the regularity arises naturally only 
when a three-body boson interaction is included. This 
result points to the origin of the three-body boson inter- 
action, suggesting the optimal IBM description of 7-soft 
nuclei. The principal results presented in this Letter do 
not depend on details of the EDF, and suggest us a com- 
prehensive picture of triaxial shapes of atomic nuclei in 
a fully microscopic way, including a solution to the long- 
standing problem of the energy-level pattern of odd-spin 
states. 
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